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Abstract 

We use continued fractions for a study of the thermodynamic properties of the periodic nonuniform spitt- 
le isotropic XY chain in a non-random/random (Lorentzian) transverse field. The obtained results permit 
to examine the influence of a magnetic field and randomness on the spin-Peierls dimerization. 
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Spin- 2 XY chains provide an excellent ground for a rigorous study of different properties of low-dimensional 
quantum magnetic systems. By means of the Jordan- Wigner transformation such spin models can be mapped 
onto noninteracting spinless fermions and many statistical mechanics calculations can be performed exactly 1 
In the present Letter we shall exploit the spin-i isotropic XY chain in a transverse field (the transverse XX 
chain) for an analysis of the effects of periodic nonuniformity on the thermodynamic properties. 

Starting from the late 60s a number of authors considered the spin-i XY chain in a transverse field on two 
sublattices. 2 ' 3 Later it was found that the XX chain (without transverse field) represents a simple system with 
a spin-Peierls instability. 4 Some further studies of the spin-Pcierls instability in the spin-i XY chains can be 
found in Rcfs. 5, 6 (within the adiabatic limit) and Ref. 7 (within the non-adiabatic limit). In what follows 
we want to emphasize that the thermodynamic properties of the periodic nonuniform transverse XX chains 
can be examined in a very general fashion including into consideration even regularly alternating randomness. 
Such an approach is based on the continued-fraction representation of the (random-averaged) diagonal one- 
fermion Green functions. The involved continued fractions can be calculated exactly for any finite period 
of regular nonuniformity that immediately yields rigorous explicit expressions for the diagonal one-fermion 
Green functions, the density of magnon states, and hence for the thermodynamic quantities. We show that the 
periodic nonuniformity may induce a number of interesting properties, e.g. plateaus in the low-temperature 
magnetization vs field dependence or a nonzero magnetization in the zero average field. As an example we 
consider in some detail the chain with the random Lorcntzian transverse field in the case of the period 2. The 
obtained density of states for such a chain permits us to extend the consideration of the influence of an external 
field on the stability of the spin-Peierls dimerized phase 5 discussing in more detail the non-random case and 
reporting for the first time the effects of randomness. The interest in the properties of the spin-Peierls systems 
considerably increased since the discovery of the inorganic spin-Peierls compound CuGeOs. 8,9 The behaviour 
of the spin-Peierls systems in the presence of an external field or randomness attracts much interest both from 
the experimental and theoretical viewpoints. 9 We believe that the discussed possibility to analyse rigorously the 
influence of non-random/random field within the frames of the simple model may be of interest for the theory 
of spin-Peierls compounds. 
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We consider N — > oo spins i on a circle with the Hamiltonian of the nonuniform transverse XX model 

N N 

H = J2 OX + 2 ^ + (1) 

n— 1 n— 1 

where the transverse fields f2 ra are assumed to be independent Lorentzian random variables each with the 
probability distribution p(Q n ) — ^ ^ -n^^+r 2 ' ^ e assume regular nonuniformity, i.e. the mean value of the 
transverse field f^on at site n, the width of its distribution r„, as well as the exchange coupling between 
neighbouring sites n and n + 1 vary regularly from site to site with a period p, i.e. the sequence of parameters 

is f^lTl/lffer^ • ■ ■ ^0p^pIp^Ql^lIl^02^2l2 ■ ■ ■ QopFpIp .... 

Our task is to examine the thermodynamic quantities of the spin model (|l|). Making use of the Jordan- 
Wigner transformation 1 one gets a tight-binding model of spinless fermions on a circle. Introducing the tem- 
perature double-time one-fermion Green functions and averaging a set of equations for the Green functions 
(because of the Lorentzian probability distribution this can be done exactly 10 ) one can further make use of 
the continued-fraction representation for the diagonal Green functions exploited earlier for similar models by a 
number of authors. 11-14 For any finite period of varying the mean value of the transverse field, the width of its 
distribution and the exchange coupling the involved continued fractions are periodic and can be evaluated by 
solving quadratic equations. Thus one gets rigorously the diagonal Green functions and therefore the density of 



magnon states p(E) (a bar denotes the random-averaged quantity), which yields the thermodynamic quantities 
for the spin model ([!]). Some examples of such calculations can be found in Ref. 15. Let us present the result 
for the density of states of a chain having period 2 that will be used afterwards 
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In the non-random case Ti = T 2 = Eq. (0) yields a nonzero value of the density of states p(E) = 
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nly for ~B'(E) = 41$ 1$ - [(E - fl 01 )(E - fi 02 ) - l\ - if] 2 > (two subband magnon en- 



ergy spectrum). In the uniform case f2oi = ^02 = ^0, Ii = T2 = T, Ji = I2 = I Eq. (g) may be rewritten as 



p(E) = + -Im . 1 . that coincides with the result reported in Ref. 16. 

Let us briefly discuss some special properties of spin model (|l]) induced by regular nonuniformity. Con- 
sider at first the non-random case T n = 0. The main result of introducing the periodic nonuniformity is a 
splitting of the initial magnon band into several subbands. The number of subbands does not exceed the 
period of the nonuniformity. The splitting of the magnon band into subbands has several remarkable con- 
sequences for the thermodynamic properties. For example, the dependence of the transverse magnetization 
00 

m z = — A f dEp(E) tanh on the transverse field £! (f2on = ^0 + ^[j n , ^o„ are fixed) at low temperatures 

—00 

is composed of sharply increasing parts separated by plateaus, the number of which is determined by the number 
of subbands. This can be nicely seen in Figs, la, lb for the case of a chain with p = 4. 

One of the interesting magnetic properties of the periodic nonuniform spin-^ transverse XX chain is a possi- 

JV 

bility of a nonzero transverse magnetization m z at zero average transverse field ^2 ^0n = 0. A similar property 

71=1 

was found for the spin-i transverse XX chain with correlated disorder. 17 To illustrate this let us consider a chain 
having the period 4 in which the site n + 1 with the transverse field fl' 02 < is surrounded by the strong couplings 
1 1 = I2, whereas the site n + 3 with the transverse field — fl' 02 > is surrounded by the weak couplings I3 = I4 = 
0. The transverse fields at site n and site n + 2 are assumed to be zero, i.e. £l' Q1 — il' 03 = 0. Naively reasoning one 
may expect that the local transverse magnetization at site n+1 has smaller value and opposite direction with 
respect to the local transverse magnetization at site n + 3 and therefore a nonzero total transverse magnetization 
at zero average transverse field may be anticipated. The explicit expression for the density of states for such a 



chain reads p(E) = i \S ( E - i(ft 02 - y/Q* 2 + 87?) J + 6 (E) + S ( E - i(ft 02 + ^02 +8/?)] + 6 (E + ft 02 ) 

JV 

(solid curve in Fig. la) and for T = one gets m z = — | ^ at ^2 ^On = (solid curves in Figs, lb, lc). The 

71=1 

picture becomes more complicated if the weak couplings I3 = I4 have small but nonzero value (see Fig. 1). In 

N 

that case for T — we have m z = at ^ Qo n = 0, however, the position of the magnon subbands provides 

71=1 

an interesting dependence of m z on temperature (dashed curve in Fig. lc). The latter dependence reminds the 
'order from disorder' phenomenon, 18 i.e. increasing of order with increasing temperature. 

The main effect of the diagonal Lorentzian disorder is a smearing out of the band structure. The details of 
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the smoothed magnon subbands may be essentially different in the case of the uniform disorder (all r„ are the 
same) and the nonuniform disorder (T„ arc different). In the latter case one may observe a different degree of 
smearing out of the different subbands that results in a different degree of smoothing of the different 'steps' in 
the step-like dependence of on f^o • 

Finally, let us discuss the spin-Peierls dimerization in the presence of non-random/random field following 
Ref. 4. For this purpose we introduce the dimerization parameter 5, < 5 < 1 assuming that = |/|(1 + 5), 
I -^2 1 = The magnetic ground state energy per site immediately follows from the derived density of states 



eo(<5) = —\ J dEp(E)\E\, the elastic energy per site is ad 2 and hence the total energy is £(5) ~ eo(S)+a8 2 . 



We shall study the change of the total energy £ (S) — £(0) as a function of the dimerization parameter S revealing 

in such a way the instability of the chain with respect to dimerization. For the non-random case eo(<5) can be 

i> 

expressed through the elliptic integral of the second kind E(?/>, a 2 ) = J dipy/l — a 2 sin 2 <p, namely, 

o 



eo{d) = {*> ( n 01 - ^ 2 ) 2 + lei 2 ) lQai + Qml U - T, 



0, if ^Shx - «o 2 ) 2 + 16/ 2 < |O 01 + fi 02 |, 



ib = < 



arcsin J M^ggjgf , if V("oi - « 02 ) 2 + 16/ 2 o 2 < |fi i + fi 02 | < y/Wi - ^) 2 + 16/ 2 , ( 3 ) 



% f, if Ifioi + «o 2 | < - O 02 ) 2 + 16/ 2 <5 2 . 

Eq. (|J) in the uniform limit fioi = ^02 = yields the result obtained in Ref. 5. We shall also consider the 
solution S* of the equation d£ SP = that may have a nonzero solution 5* 0. In the non-random case the 
latter solution may arise as a solution of the equation 



V(0 01 -Oo 2 ) 2 + 16J 2 / / 16J 2 (1-Q \ / JWM_, . 
4tt(1 -o 2 ) ^ ^' (ft i ~ tt 02 ) 2 + RSI 2 / ^' (ft 01 - O 02 ) 2 + 16/ 2 H 

where F(-0,a 2 ) = f dip , == is the elliptic integral of the first kind. 
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Let us discuss the dependence £(5) — £(0) vs <5 which manifests the spin Pcierls dimerization. We restrict 
our discussion to the uniform case f2oi = ^02 = ^0, Ti = T 2 = T. Consider at first the non-random case. Note 
that for strong transverse fields |f2o| > 2|J| from Eq. (||) it immediately follows that £{5) = — ||f2o| + c-5 2 and 
hence £(S) — £(0) exhibits a minimum only at S* = 0. This tells us that in strong enough field the uniformed 
state is energetically favoured over the dimerized one. Consider further weaker fields |fio| < 2|/|. From Eq. (||) 



it follows that £(8) is influenced by the field for < 8 < ^rjj- but £(8) does not feel the presence of the field for 
]5yl < 8 < 1. On the other hand one can calculate the r.h.s. of Eq. (Q) varying 8 from to 1 finding in such 
a way the values of the dimerization parameter at which the total energy for a certain lattice (characterized by 
a) exhibits an extremum. As it follows from Eq. (^) ^\ \ as 8 ~~ * 1 that means that for < \ the total 
energy £(8) does not have the extremum of interest for < 8 < 1. Hence, we shall exclude the soft lattices 
having small a from further consideration. The dependence £{8) — £(0) vs 8 at different f2o is shown in Fig. 2a. 
One finds that at fio = the quantity £ (6) — £(0) exhibits a minimum at 5* ^ 0. The position of this minimum 
does not change switching on the field. However, after a certain value of the transverse field (for which Eq. ([|) 
has an additional solution 6* = 0) an additional local minimum in £($) — £(0) appears at 8* = 0. The two 
minima are separated by a maximum at an intermediate value of the dimerization parameter. With further 
increasing of f^o the depths of the minima become the same at a certain value of the transverse field and then 
the minimum at 8* = becomes the global one. Finally at a certain field (for which Eq. (|J) has the solution 
8* = Tpjy) the minimum at the nonzero value of the dimerization parameter abruptly disappears that means a 
complete suppression of the dimerization by the field. The above sketched picture represents a typical scenario 
of a first order phase transition and is illustrated by the phase diagram in the f2o — ol plane shown in Fig. 3a. 

The dependences of £{S) — £(0) vs 8 in the presence of randomness at f^o = and ilo 7^ are shown in Figs. 
2b, 2c. In the case f2o = with increasing of T the position and the depth of the minimum at 8* 7^ continuosly 
decrease and for a certain value of the strength of disorder the minimum occurs already at 8* = 0. That 
mean, that a random field with zero mean value can suppress the dimerization, too. However, the dimerization 
parameter 8* vanishes according to a second order phase transition scenario in contrast to the transition in a 
non-random field. The influence of the randomness on the dependence £(8) —£(0) vs 8 at Slo 7^ can be traced 
in Fig. 2c. The phase diagrams in the plane T — a are shown in Figs. 3b, 3c. 

It is necessary to mention that we have considered here only the stability of the dimerized phase against 
the uniformed one. Though, the qualitative picture of a first order phase transition in a uniform field may 
be correct, we notice that both experiments and approximate analytical treatments as well as exact numerical 
computations (mainly for CuGeOa) show a transition from the dimerized to incommensurate phase. 9 Evidently, 
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assuming the simple ansatz for the lattice distortion 61626162 ■ ■ ■ , 61 + 62 = we were able to compare the 
ground state energies only for the dimerized and uniformed phases. To detect a transition from the dimerized 
phase to the incommensurate phase with increasing of the field one must analyze the ground state energy of a 
chain having sufficiently large period (say, a chain with p = 12, for which one can still easily obtain the density of 
states 19 ). In the presence of random fields (or/and random exchange couplings) even more complicated lattice 
distortions should be examined and the elaborated approach for the calculation of the ground state energy of 
the nonuniform spin-i transverse XX chain provides some possibilities to perform such an analysis. 

To summarize, we have suggested a general scheme for the study of thermodynamics of the regularly al- 
ternating spin-i transverse XX chains. It is based on the continued-fraction representation for the diagonal 
one-fermion Green functions. We have discussed briefly some magnetic properties of the periodic nonuniform 
spin-i transverse XX chains such as plateaus in magnetization curves or a nonzero magnetization at the zero 
average field, as well as the influence of (Lorentzian) randomness in the transverse field on these properties. We 
have examined also the stability of spin-Peierls dimerized phase in the spin-i XX chain in the presence of non- 
random/random (Lorentzian) transverse field using the rigorous expression for the (averaged over randomness) 
ground state energy of the corresponding alternating chain. 

The present study was partly supported by the DFG (projects 436 UKR 17/20/98 and Ri 615/6-1). It 
was presented at the International Workshop and Seminar on Cooperative Phenomena in Statistical Physics: 
Theory and Applications (Dresden, 1999). O. D. is grateful to the Max-Planck-Institut fur Physik komplexer 
Systeme for the hospitality. He is also indebted to Mrs. Olga Syska for continuous financial support. 
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FIGURE 1. Magnetic properties of the chain having period 4. The local transverse fields are fion = flo + ^'om 
Q' 01 = ^q 3 = 0, n' Q2 = —^04 = — 1) the exchange couplings are |7i| = |/ 2 | = 0.5, I-Z3I = I/4I = (solid curves), 
I/3I = I/4I = 0.05 (dashed curves), I/3I = I/4I = 0.25 (dashed-dotted curves), I/3I = I/4I = 0.5 (dotted curves), 
a) Density of states, b) Transverse magnetization vs transverse field fi at T — 0. c) Transverse magnetization 
vs temperature at Qq = 0. 



FIGURE 2. £{8) - £(0) vs S for a chain with |/| = 0.5, a = 0.5. a) Non-random case r = 0, fi = 
0, 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, 0.07, 0.08, 0.09, 0.1 (from bottom to top), b) £1 = 0, T = 0, 
0.001, 0.003, 0.006, 0.01, 0.02, 0.03, 0.05 (from bottom to top), c) fl = 0.04, T = 0, 0.001, 0.003, 0.006, 
0.01, 0.02, 0.03, 0.05 (from bottom to top). 

FIGURE 3. Phase diagram of the spin-i transverse XX chain with |/| = 0.5 as it follows from the 
dimerization ansatz for the ground state energy, a) Non-random case T = 0. A — dimerized state. Bi, B 2 - 
both dimerized and uniformed states are possible, moreover, in Bi the former one is favourable, whereas in B 2 
the latter one is favourable; on the line that separates Bi and B 2 the depth of both minima in £{5) — £(0) is 
equal. C — uniformed state, b) Random case with mean value of the field Oo = 0. A — dimerized state. C — 
uniformed state, c) Random case with mean value of the field f^o = 0.04. A, Bi, B 2 , C arc of the same meaning 
as in a). 



9 




10 



0.001 



0.000 



as)^(O) 

0.0000 
-0.0002 

0.0001 
0.0000 



0.00 



